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The self-consistent propagation of generalized Di [coherent-product] states and of a class of 
gaussian density matrix generalizations is examined, at both zero and finite-temperature, for arbi- 
trary interactions between the localized lattice (electronic or vibronic) excitations and the phonon 
modes. It is shown that in all legitimate cases, the evolution of Di states reduces to the disentan- 
gled evolution of the component D2 states. The self-consistency conditions for the latter amount to 
conditions for decoherence-free propagation, which complement the D2 Davydov soliton equations 
in such a way as to lift the nonlinearity of the evolution for the on-site degrees of freedom. Although 
it cannot support Davydov solitons, the coherent-product ansatz does provide a wide class of exact 
density-matrix solutions for the joint evolution of the lattice and phonon bath in compatible sys- 
, terns. Included are solutions for initial states given as a product of a [largely arbitrary] lattice state 

and a thermal equilibrium state of the phonons. It is also shown that external pumping can produce 
' self-consistent Frohlich-like effects. A few sample cases of coherent, albeit not solitonic, propagation 

\wJ , are briefly discussed. 
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Davydov's model of soliton propagation in molecular chains is not rigorous but its transparency and physical 
appeal continue to encourage intense work regarding its utility as a practical approximation, at both zero and finite- 
temperature ||. Alternatively, it is commonly hoped that more elaborate Hamiltonians and/or wave functions may 
render the model exact or improve its accuracy, and various proposed refinements have attempted to do so Q] . From 
the latter point of view, it is quite intriguing that the conventional approach to the problem stops short of probing 
the self-consistency conditions for the wave function ansatz. Indeed, the Davydov model relies on two fundamental 
assumptions: i) A D± [D2] ansatz state is a solution of the Schroedinger equation for the Davydov Hamiltonian. ii) 
The time-dependence of the ansatz state can be obtained by treating the ansatz parameters as canonical variables in 
a Hamilton's functional given by the average of the Davydov Hamiltonian on the ansatz state. The second conjecture 
finds a self-consistent foundation in the variational principle of least action (5), while the first is suspended, for 
C7^ practical purposes, and the state ansatz is cast as a variational trial ansatz. But the ansatz state is likely to be a 
good [dynamical] trial state in those situationsthat are 'close', in some suitable, perturbative sense, to a self-consistent 
model for which this state is an exact solution. Therefore, one can start, conceivably, by questioning what particular 
circumstances allow an entangled superposition of coherent or thermal-coherent [Gaussian] phonon states to preserve 
its cohesion in time. 

The present paper aims to address this problem in a general setting, at both zero and finite-temperature. The 
outcome, which can be viewed, eventually, as a generalization of two previous theorems, due to Brown 0, on the 
validity of the standard D2 and D\ states, has the unexpected effect to place the Davydov ansatz problem in a 
fresh perspective. Our starting point is a generic, not necessarily unidimensional, lattice of N interacting monomers, 
each coupled in turn to a common phonon [bosonic] band. No a priori assumptions are made as to the nature 
of the lattice [monomer] degrees of freedom [excitonic or vibronic], of the phonon modes [acoustic/longitudinal, 
optical/dispcrsionless, etc.], or of the site-to-site and lattice-phonons interactions. Not surprisingly, however, the very 
nature of the coherent states, as states specific to harmonic systems, limits the type of appropriate lattice-phonon 
interactions to bilinear terms in the phonon degrees of freedom at zero temperature, and to only linear terms at finite- 
temperature. This notwithstanding, the critical constraints concern the accompanying, entangled lattice states and 
the lattice contribution to the lattice-phonon interaction. Somewhat contrary to the widely held view that the coherent 
phonons should drive the (self-trapped) lattice configuration, it turns out that phonon coherence is also essentially 
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conditioned by the nature of the entangled lattice states. We find that self-consistency requires, both at zero and at 
finite-temperature, that the dynamics of D\ superpositions be reduced, remarkably, to the disentangled propagation 
of the Z?2 components. Therefore no self-consistent D\ model can generate soliton equations coupling the separate Di 
states. Further, the lattice state in any D2 product must satisfy constraints leading to [slightly modified] decoherence- 
free evolution As a result, the self-consistent lattice states propagate unitarily and all potential nonlinearities 
in their effective equation of motion cancel identically. Hence, no D2 model based on a linear decomposition of the 
lattice state can generate Davydov-type soliton equations for the associated [linearly independent] lattice amplitudes. 

A stronger and quite peculiar result concerns systems with lattice-phonon couplings linear in the phonon coordinates. 
First, only such systems can support generalized, Gaussian D2 states at finite-temperatures. Second, the corresponding 
self-consistency conditions for lattice states become precisely the constraints that typically define decoherence-free 
subspaces [DFS], as recently considered in connection with noiseless quantum computation ||, [Q, and the effective 
Hamiltonian reduces to the unperturbed lattice Hamiltonian [up to a phonon- modulated energy shift, which can be 
transferred onto the phonon state]. When the coupling is time independent, the associated phonon displacements 
perform simple harmonic oscillations around a displaced equilibrium position. One is so led to the conclusion that 
finite-temperature coherent-product propagation is possible if and only if the lattice-phonon coupling is linear in the 
phonon coordinates and decoherence-free subspaces exist for the lattice subsystem. 

From the point of view of the soliton problem, this outcome relates to a rather self-evident idea in the context 
of decoherence-free propagation. That is, if an unperturbed lattice can support soliton states that evolve entirely 
within a DFS, such solitons will propagate unperturbed [at arbitrary temperatures], regardless of the state of the 
environment. Obviously, coherent-product states loose any special significance in such a solution, and the problem 
falls outside the scope of the present paper. 

Nevertheless, the self-consistent coherent-product ansatz demonstrates a sufficient number of notable features to 
remain attractive by itself. For instance, let us recall that null phonon displacements place the bath in a thermal 
state, hence initial null displacements reduce any initial Gaussian D\ state to the ubiquitous product of a lattice 
state and a phonon thermal state. Since the initial lattice state can be an arbitrary distribution on a direct sum of 
[orthogonal] DFS, we obtain the following corollary, with reference to the theory of decoherence-free propagation. 

Let a system interact with a boson bath through a coupling linear in boson coordinates. Any distribution on a direct 
sura of system DFS develops, when brought in contact with the bath in a thermal state, into a [strongly entangled] 
Gaussian D\ state. 

Because the expression of the D\ state is exactly known, the development of the system-bath entanglement in 
such a model can be monitored precisely. The result is a remarkable counterexample to the standard picture of 
relaxation through thermal contact. Regardless of any intrinsic characteristics, the bath is forced into a coherent, 
nonequilibrium state, while the DFS components of the system state evolve in an unperturbed manner. This confirms 
earlier claims that decoherence may produce in fact coherent states H, and therefore caution should be exerted in 
assuming that a heat bath remains at all times in thermal equilibrium |9| . It also corroborates the similar conclusion 
of a recent study |l(J of decoherence in a quantum system interacting with a [macroscopic] measuring apparatus. 
When the correlations between distinct DFS vanish in time, the system is driven toward a statistical superposition of 
decoherence-free, pointer states and we retrieve a typical example of environment-induced superselection Jjyj]. But in 
another interesting limit, which arises under time independent interactions, the bath modes can be prepared such that 
the initial distortions match the displaced equilibrium positions for the coherent oscillations. In this situation, the 
bath remains in a stationary nonequilibrium state, while the system evolves unitarily, according to the unperturbed 
dynamics [up to time independent DFS energy shifts], irrespective of any entanglement with the bath. Leaving aside 
questions of stability, such a state provides an apparent counterexample to environment-induced superselection, within 
the same physical system. 

Along a different line of inquiry, proper external pumping can be used to promote decoherence-free propagation 
in systems that otherwise may not have the necessary symmetries. We find, incidentally, that such a process can be 
accompanied by a Frohlich-like effect on the bath modes. That is, pumping at a frequency attuned to the lattice 
subsystem may result in a macroscopic displacement of a bath mode of a different frequency, while other reservoir 
modes remain in thermal equilibrium. 

The paper is organized as follows. We begin with the familiar pure-state case in Sec. II. Section IV develops the 
density-matrix generalization, which includes Davydov's thermal ansatz and allows a straightforward approach to the 
finite-temperature problem. To this end, we employ a formal framework, outlined in Sec. Ill, based on the 'square- 
root' decomposition of the density-matrix and the concepts of thermal vacuum and thermal Fock space introduced 
in thermofield dynamics [TFD] (l2| . This formalism can be regarded in effect as a version of TFD without auxiliary 
systems. Our choice is motivated by the notable technical advantage that pure state calculations can be effortlessly 
transcribed into the density-matrix domain, with a minimum of adjustments. Section V examines and discusses 
some popular versions of Davydov's model, alongside with sample self-consistent examples, including a case with 
time-dependent interaction (external pumping). A summary and concluding remarks are provided in Sec. VI. 
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II. SELF-CONSISTENT DYNAMICS OF Di STATES AT T = 

Let the lattice-phonon Hamiltonian be 

H = H lat + H° ph + W , 



where 



N-l 



Hlat = E £ n C ti C n + Via 



n=0 



describes the iV-site lattice, with Vi a t subsuming all hoping interactions, 

H° h = J2^lb q 



(1) 



(2) 



(3) 



is the bare phonon Hamiltonian, and W accounts for any lattice-phonon and/or anharmonic phonon interactions. Let 
us search for self-sustained soliton states |^) in the slightly generalized Davydov D\ ansatz 



i*> = X>>i/j«> , 



(4) 



where the |a)'s are orthogonal, but not normalized, lattice (vibronic) states, (a'\ a) ~ 5 a i a , ^2 {a \ a) = 1, and the 
|/3q,)'s are coherent phonon states, 



exp 



|o) 



ph 



(5) 



Here the |a)'s are allowed to contain local excitations in different numbers and are not necessarily confined to the 
same site. The [9>) ansatz above is, thus, general enough to cover both single-quantum and multiquanta solitons, as 
well as their eventual superpositions. The set of occupied states {|a)}( Q |^^ is to be regarded as embedded into an 
orthogonal basis {\a)} of lattice states; the unoccupied \a)'s, (a | ^) = 0, may be assigned by default |/3 Q ) = |0) h . 
For each displaced phonon vacuum \(3 a ), let us construct the corresponding orthonormal Fock basis 



n4 



n q \ 



Pqa) \fia 



(6) 



where |. . .n q>a ■ ■ ■) carries n q . a displaced quanta of mode q relative to \f3 a ). In particular, \l q . a ) — (b\ — qa ) \/3 a ) 
denotes the first excited state of the displaced mode q relative to the displaced vacuum \j3 a ). The set 
{\a) \. . . n q ^ a •••)}<, Y]n >o obviously provides an orthogonal basis for the overall lattice-phonon system. 

q 

Extending the idea originally applied by Brown in Ref. M, the conditions under which the state |\&) of Eq. (Q) 
is compatible with the dynamics driven by H will be derived from the expansion of the corresponding Schroedinger 
equation 



(J) 



in this displaced phonon basis. Indeed, substitution of the explicit form of followed by some straightforward 
algebraic manipulation, leads to 



E l d > + \ E " q J* qa ) \a) \0 a ) + E 4o 



la) 1 



a \ q / a, q 



(8) 
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Since the time derivative of |^) on the left-hand side of Eq. fl|) carries only terms in displaced vacuums and their 
first excited states, it is foreseeable that the main constraints will emerge from the requirement that higher order 
contributions on the right-hand side vanish identically. As previously hinted, we can expect the interaction to 
be limited to bilinear terms in the phonon coordinates. But it will also become apparent that the form of the 
effective equation of motion and the self-consistency conditions for the lattice states |a) are shaped by this same 
requirement. So consider first the contraction of both sides of Eq. (||) with states \a') \. . .n q ^ a i . . .) carrying more 
than one phonon quantum, ^2n q > 1. After conveniently expanding the contribution in W along displaced phonon 

q 

bases and rearranging the terms, we are left with 

[(a' \H ex \ a) + (a' \(J3 a \W\ f3 a ) \a)-ih (a' \ a)} (. . . ry , Q , . . . | /3 a )+ 



(«' |(l,,a|W| | a) {...n q >, a > ... | 1 



q.Oij 



+ J2 J2 {a' \{...m q}a ...\W\/3 a )\a){...n q ', a >...\...m q}Ct ...} = 0. (9) 

a {m q } 
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Under the reasonable assumption that the form of the interaction W limits the last sum above to a finite number of 
terms, a nontrivial solution will be compatible with the infinite number of constraints (|^) if and only if all quantities 
multiplying the nonvanishing overlap factors (. . . n q i, a i ... | ... k q ^ a . . .) cancel identically. In particular, the terms 
under the last sum in Eq. (^) imply 

(...m q>a ...\W\ p a ) =0, J2 rn «'<* >h ( 10 ) 
Hence W can only contain terms in b q , b q and b q b q i , and must be of the form 

w=j2(K w i+ b i w t)+J2 b l b i' u ^ > ( n ) 

9 q,q' 

where the operators w q and u qq i act on on-site degrees of freedom only, as (polynomial) functions of c„, c+, and the 
hermiticity of W requires u qlq = u qq / . Further, from the second sum in Eq. (|^) it can be inferred that 

(a' | (l g , Q I W\ a ) | a) = V«, a ^ a' . (12) 

But since |a') spans an entire lattice basis, the above constraints imply in fact that the occupied states \a) must 
satisfy eigenvalue equations of the form 

G q , a | a) = 7,, Q \ct) , Va, q , (13) 
for G qM = (l q: a \W\ f3 a ). In the ansatz ( [Tl| ) for W, the operators G q _ a read 

G«j,q = w q + ^ I3 q > a u qql . (14) 
q' 

Similarly, from the first sum in Eq. (|^) it follows that 

{a'\Hj\ a) + (a' \(/3 a \W\ (3 a )\ a) -ih(a' \ a) = , Va' ^ a , (15) 
which shows that the occupied |a)'s, (a \ ^ 0, should evolve according to an effective Shroedinger equation 

ih \a) = [H ex + (J3 a \W\ /3 a ) + fl a ] \a) , (16) 

with Q a a scalar functional, and 
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(0 a \W\0 a ) 



E 

Q 



(0* qa w q + qa w\) + 0* qa q'a u qq' 
q,q' 



(17) 



It is already evident, from Eqs.(13) and (filf), that the evolution of distinct lattice states \a) is reciprocally decoupled, 
unless the eigenvalues fj, qja exhibit a dependence on some qa ' with a' ^ a. It will be seen shortly that this is not the 
case. At the same time, it can be recognized that the eigenvalue equations ( |l3| ) implement an effective decoupling of 
the lattice-phonon interaction, very much in the manner of the effective decoupling responsible for decoherence-free 
subspaces ||. The latter process requires that, for a lattice-phonons interaction of the general form ^ W^ at ^ , 

k 

the lattice state \^i a t) be such that wj. lat ^ \^i a t) — Mfe \^lat) for all k, and at all times. A subspace of lattice states 
complying with these constraints is termed a decoherence-free subspace (DFS). Obviously, states belonging to a 
lattice DFS are completely decoupled from the phonon dynamics and evolve unperturbed. A closer examination of 
our conditions ( |l3| ) shows that, in fact, the present process differs from a true DFS selection simply by the assumption 
of a coherent ansatz for the state of the phonon modes. When the interaction in Eq. ( |Tl| ) is applied to any product 
| a) \0 a ), the result reads 



W \a) \0 a ) 



^ Pq'aUqq' 

q' 





\p*) + 


**TPq a wl \a) 








- q 





and it is immediate that condition (13) is just the proper condition for the decoupling of the first term on the right- 
hand side. It must be kept in mind, however, that we have arrived at this result without any assumptions on the form 
of W or on the separability of the D2 terms at the outset. 

The details in the equation of motion (|l^) and the self-consistency constraints (|l3|) will follow, as can be anticipated, 
from a balance of the left-hand side of Eq. (|J) with similar, nonvanishing terms on its right-hand side. Indeed, let us 
first contract Eq. (ph with states \a') \0 a ')- Taking into account conditions (Eq) above, we are led to 



ih (a | d> = (a \{H ex + (0 a \W\ a ))\ a) + £ Uu q \0 qa \ 2 - % - (p qa qa - q J* qa ) 

q 

which, in conjunction with Eqs.([l^), identifies 

= E 1^9" | 2 ~~ ~2 {P^Pqa ~~ ftqaPqaj 

q L 

Likewise, contracting Eq. (@) with states \a') \l q , a ') yields 

ifipq'a' - Tiu) q >P q > a > {a' I a')<5 ,(Q'|*) = ( a ' l(V,a' \W\ a >)\ a') 5 o, (a'|*) + 



a a) 



(18) 



(19) 



+ {(a'\H ex \a) {l q ,, al I a ) + (a' |(1,,, Q , \W\ a )\ a) - ih (a' | d) (l q ', a > | a )} . 

Using again ansatz (^Tj) for the interaction W and Eq. ( p"2| ) in the form (a 1 \G q)OL \ a) ~ <5 a / Q , we note that 
]T (<*' KV,a' \W\ a )\ a) = Y [< 1 «'.a' H\ M («' |G,, a | a) + qa {a' \w\\ a) (l q ,, a , \ a )] 



(20) 



Y P qa (<*' K| a) (l q ,, a , \0 a )=Y (<*' \{Pa \W\ a )\ a) (l q ,, a , \ a ) 

ol^ol' ,q a^a' 



(21) 



It becomes apparent now that the sum on the right-hand side of Eq. (|20j) cancels by virtue of conditions (|l^) , while 
the remaining expression complements Eq. (113) , and shows that the exact form of Eqs.(Q3h is 
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Gq. a | a) = ih$ qa - huj q l3 qa I a) , Va, g 



(22) 



Finally, it must be verified that the occupied lattice states remain orthogonal at all times [(a (t) | a' (t)) ~ <5 Q ,a'] 
under the propagation described by Eqs.(|lq). Indeed, it is immediate that for any two occupied states, (a | \&) 7^ 0, 
<a'|*)#0, 

ih [(a' I d) + (a' I a)} = (a'\ (p a \W\ f3 a ) \a) - (a'| (/3^ \W\ M |a) , (23) 
which, in view of Eqs.(|l7|) , ( |22| ) and (fbt), becomes succesively 



51 [/W«'<*' (a' a) - &* qa ,fi q > a (a' \n qq >\ a)] = ^T^faa 



a) - (a' |it M '| a) 



(24) 



The last line above follows from the condition that u qlq — u qq > in order to secure the hermiticity of W. 



To sum our results to this point, a (generalized) D\ state |W) [Eq. 
dynamics if and only if 



describes a self-consistent lattice-phonon 



a) the interaction W is of the form (11) 



b) the lattice states \a) are driven by the effective Schroedinger equations (jig ) and each satisfy, simulteneously, a 
(large) number of DFS-like constraints imposed by the eigenvalue equations (|22| ) for the operators Gq, Q ; 

c) the phonon displacement parameters (5 qa obey the evolution equations 



(25) 



where 7 g]Q , denotes an eigenvalue of G g , a (assuming any exists) and we have indicated explicitly the dependence on 
the set of (3 qa 's and on the interaction factors w q and u qq i . 

As suggested earlier, self-consistency is seen to require that the propagation of individual |a)'s and their associated 
phonon parameters {(3 qa } be decoupled from similar terms. This rather strong result shows that, under quite general 
conditions, a lattice-phonon system can support D\ states if and only if it can support the component orthogonal D2 
states I^q) = \o) \Pa) individually. In particular, it ensues that a given system can sustain standard D\ states 



N 



l*l> 



H=nC\ \(3 n ) |0) e 



(26) 



if and only if it can sustain on-site, single-quantum coherent states = |0) |/3 n )- Even under such circum- 
stances, it turns out that the D\ states can only yield static lattice distributions [static 'solitons']. Indeed, according 
to Eq. d2~i|), (a I a) — const [consistent with the unitary evolution of \a)] and, correspondingly, the probabilities 
\tp n \ 2 in any self-consistent, standard D\ state [Eq. (^6|)1 are conserved in time. The same holds true for higher 
order (multiquanta) D\ states based on on-site |a) states, as well as for their superpositions. As a corollary, mobile 
lattice distributions can arise if and only if the driving Hamiltonian accomodates D2 states with delocalized lattice 
components. 

Let us now note that, while Eqs.(|l6|) are just the usual Davydov equations for the [D2] lattice factor states, Eqs.(|25|) 
do not involve the |a)'s themselves, hence no parameters defining these states can appear in the evolution equations 
for the f3 qa 's. In other words, a solution for the f3 qa 's, when it exists, will not show a functional dependence on the 
parameters of \a), and substitution of the explicit expressions for the /3 ga 's into the Hamiltonian of Eq. ( |l9| ) cannot 
lead to nonlinear equations for |a). 

To clarify the relation between this result and Davydov's soliton equations, rewrite the equations of motion for the 
P qa 's as [Eq. (H), H a) = 1] 



ih/3 qa - huj q l3 qa = (a\ G 9iQ \a) . (27) 

A linear parametrization of \a) in terms of some fixed lattice basis states, \a) = X] a fc w ^ u induce, apparently, 

k 

nonlinear equations typical of Davydov's theory. In fact, it can be ver ified that Eqs.© and @ do recover the 



G 



£>2 form of the soliton equations for the appropriate Hamiltonian and choice of |a). However, Eqs.(|l6|) and rtgj] ) 
cannot describe a self-consistent dynamics unless complemented by the requirement that |a) also be an eigenfunction 
of Gq t0l . But, in that case, the average on the right-hand side of Eq. ( p7| ) reduces to the corresponding eigenvalue, 
and the propagation of |a) becomes linear. It is interesting that although the exact form of Eqs. ( |22|) is known in the 
literature for certain extensions of Davydov's model [see, e.g., Ref. pfl ], their interpretation was merely as a particular 
type of nonlinear equations for the phonon parameters. Of course, the linear effective evolution of the lattice states 
|a) does not preclude, by itself, nontrivial phenomena in the lattice dynamics. Provided appropriate solutions exist, 
it is conceivable that nonlinear effects can yet arise, and be modulated by the phonon contribution. However, the 
circumstances leading to such effects are, obviously, rather restrictive. 

The constraint that \a) be a (simultaneous) eigenfunction of the operators G qa can be understood also as a 
restriction on the allowed initial state |a), as seen by substituting \a (t)) — U a (t) |a(0)), with U a (t) the unitary 
evolution operator corresponding to the effective Hamiltonian of Eq. (161). Further, for given a, the set of constraints 
( p2| ) is equivalent to the single initial-state constraint 

A a |a(0)) = 1 (28) 

for 

oo 

A a = [dtJ2 [Gt, a (*) - lla (*)] [G 9 ,a (*) " 7 9 ,a (t)] , (29) 
« 

with 

G q , a (t) = Ut (t) Gq t a (t) U a (t) . (30) 

In this point of view, the existence of self-consistent D± states is conditioned by the existence of an appropriate subspace 
of acceptable initial lattice states, determined by Eqs.(p8|) -(|30|). When such a subspace exists, the propagation of any 
Di state [Eq. (Q)] with initial lattice states belonging to this subspace may be described by the generalized Davydov 
equations ( |l6| ) and ( p7[ ) , b ut, as noted above, the apparent nonlinearity of the underlying lattice dynamics cancels 
identically. Equations ( |27| ) reduce in fact to Eqs.(^5(). In Sec. IV it will be shown that self-consistent propagation 
at finite-temperature further restricts the allowable forms of lattice-phonon interactions, and sets final bounds on the 
subspace of permitted initial states. 



III. DENSITY MATRIX FORMAL FRAMEWORK 



The approach of the previous section will be extended to the finite-temperature domain in a framework that 
borrows heavily from the thermofield dynamics (TFD) construct |l2| [a somewhat related formalism is detailed in 
[]13j| ]. Specifically, the total density-matrix for the lattice-phonon system is written as p — 77^, where the nonhermitian 
state operator 7 is defined up to a unitary gauge transformation of the kind 7 — > jU (t), UW — Wll = I. In other 
words, any two distinct state operators 7,7' corresponding to the same density matrix are necessarily related by 
7' = (t) . From the evolved density-matrix, written in the form 



p (t) = exp 



7 (0) U (t) ■ (t) 7 f (0) exp 



one can infer that 



7 (t) = exp 



1(0) U(t) 



satisfies the von Neumann-like equation 



c?7 

in— = H ■ 7 — 7 • K 
dt 



(31) 



with K the [arbitrary] hermitian generator of the unitary gauge factor U (t). In superoperator notation, Eq. (|3l| ) 
reads 



7 



^ = (H-K) 7 , (32) 



H,K 



= 0. 



where H is the superoperator defined by the Hamiltonian H, H7 = H ■ 7, and K is given by K 7 = 7 • K, 
As in TFD, the tilde notation is used here for the tilde conjugate of a superoperator A, linear or not, introduced as 


[A(a)] f = A(d f ) . (33) 

If A (&) = A • a for some linear operator A, Eq. j3^ ) yields A (d) = a ■ A\ which reduces to A (d) — a ■ A when A 
is self-adjoint. It is immediate that the tilde operation is distributive against the usual addition and multiplication 
of linear operators, but antilinear with respect to the multiplication by scalars. Also, tilde-symmetric superoperators 
map self-adjoint operators onto self-adjoint operators. 

Obviously, from Eq. ( |3l| ) one can retrieve the von Neumann equation for the density-matrix, which in superoperator 

form reads ifidp/dt = (tl — flj p. The gauge generator K does not bear any physical significance, but, as will be seen 

in the next section, may prove instrumental in streamlining the calculation, although it can suffer restrictions under 
a given ansatz for j. Note that while H can contain only direct superoperators generated by the usual observables 
[or creation and annihilation operators], K can only contain the corresponding tilde conjugates. 
If the space of linear operators is endowed with the Hilbert space structure induced by the inner product 

(a\6j) = Tr , (34) 

7 spans the sphere of unit norm operators [(7I7) = 1] , and the average of an observable becomes 

Tr (Op) = Tr (0*frf) = (7 |0| 7) . (35) 

An orthonormal operator basis {cr}, (<3i<r') = S aa r , determines a basis representation of 7 as 

7 = X>fiO*' (36) 

a 

and similarly for observables [superoperators]. Special attention receives the operator equivalent of the Hilbert space 
Fock basis, which implements the operator Fock space. If a and o,t are boson annihilation and creation operators 
for the Hilbert space vacuum |0), the annihilation and creation [superoperator] counterparts for the operator vacuum 
[projector] |0) (0| are a, a and, respectively, at, at. The corresponding operator Fock basis follows simply as the 
exterior product of vectors of the Hilbert space Fock basis, 

I n,m) = (at)" ( a t) m |0) <0| = _L= (at)" |0) <0| a™ = \n) (m\ . (37) 

Applying now the central idea of TFD, note that a boson thermal state operator 

7T = Z~ 1/2 cxp [-{huj/2k B T)(Ja} 
is related to the zero-temperature vacuum by the unitary, tilde-symmetric (super)transformation [Ol 

7 t = U(T)|0)(0| , (38) 

U(T) = exp [(at at - aa) 6) , 
with tanh# = exp [— huj/2kBT]. Hence, it is itself vacuum for the thermal annihilation (super)operators 

A = U (T) aUt (T) = a cosh - a f sinh 9 , (39a) 

A = U (T) aU f (T) = -a f sinh 6 + a cosh 9 . (39b) 

Just as in the zero-temperature case, the corresponding operator Fock space can be constructed on the 'thermal 
vacuum' jt by using the thermal creation (super)operators At and At. 

The advantage of using the state operator 7 and Eq. (2.2) over the density-matrix and the von Neuman equation 
comes from the evident analogy between this framework and the usual Hilbert space formalism. In addition, for the 
problem at hand, our Gaussian density-matrix ansatz (see below) has a simpler, if somewhat abstract, expression for 
7- 
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IV. SELF-CONSISTENT DENSITY-MATRIX COHERENT-PRODUCT [Di] STATES 
The zero-temperature D\ ansatz for pure wave functions (state vectors) is generalized here to the product ansatz 

7 = & ■ % ■ U (40) 



for the state operator. In direct analogy to the pure state ansatz ([g), the lattice state operators a are assumed 
orthogonal, but not normalized, in the sense of the operator scalar product (& | a') — Tr ya*a') ~ S a > a , while the 
phonon state operators ^ a are given the thermal Gaussian [coherent] form 



(41) 



with 



7g,T 



exp 



hu>, 



2Jfe B T 9 * 



(42) 



The overall conservation of probability requires ^2(a \ a) — 1. Also as in the pure state case, the a's will be considered 

a 

a subset of an orthogonal lattice operator basis {a} and the noncontributing a's, (a | 7) = 0, will be assigned by 
default null phonon displacements, (3 qa = 0. 

The ansatz (Hfl) corresponds to the density-matrix 



(43) 



which can be understood, in general, as a strongly entangled [incoherent] superposition of phonon Gaussian states. 
Indeed, averaging over the lattice degrees of freedom yields the phonon density-matrix as 



pph = Y II y exp 



TlLJn 



k B T v " 



(44) 



When the lattice factors a correspond to mutually orthogonal pure states, i.e., a — \a) (a\ U for some common unitary 
factor U, (a \ a') ~ 5 aa ', one has a'c^ = S a / a \a) (a\ and Eq. (E3h acquires the simpler and more transparent form 



P = Y ^ II J~ exp [Pi* h \ ~ Pgah] exp 



huj n 



exp [-P q J>\ + (3* qa b q ] (a\ 



(45) 



If the [unnormalized] states \a) are set proportional to the on-site, one-quantum states |1„) = cjj |0), expression (Q) 
recovers Davydov's original finite-temperature ansatz [[TBI . At zero temperature the phonon contribution reduces to 
7a — ► \Pa) (0| pft and the ansatz ( ^p| ) yields 

7 = £a-|/3 Q ) (0\ ph U . 



ph 



For lattice factors a — \a) (\\, with a common (\\, (x \ X) — 1; the state operator becomes 7 = X) l a ) \Pa) (x\ (0| 

L & 

and corresponds to a pure Dl state of the type (Q). Even more remarkably, when all displacement parameters vanish, 
Pqa = 0, all 7 Q reduce to the thermal state operator jt = Y\lq,Ti an d the density-matrix ( f43| ) acquires the familiar 

product form 



Plat ■ Pph.T 



where pi a 



£< 



and /5 p /,.,t = II (%q) 1 ex P [ — (Sw g /^T)iij!) g ] . Since the lattice state operator 



(46) 



is quite arbitrary, the ansatz (J40|) is seen to cover all situations where a product state as in Eq. (|46|) evolves into an 
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entangled mixture described by Eq. (43). Note that the ansatz ( pip] ) accounts for a much wider class of evolutions, 
because, in general, the displacement parameters (3 qa need not vanish simultaneously at any time. 



Let us now introduce, for each phonon mode q, the thermal annihilation (super)operators 

B„ = b„ cosh 9 q — h\ sinh 9„ , 



(47a) 



-b q sinh 9 q + b q cosh 9 q , 



(47b) 



with tanh#g = exp [— Tiu) q /2kBT]. Upon substituting the reciprocal transformations into 7 9jQ 
exp [/3 ?Q bJ — /3q a b g ] %,t, it is seen that j a — Y[%,a represents a displaced vacuum for B g and B g , since 



B 9 7 Q = /3 qa cosh9 q 'y a 



(48a) 



B 9 7 a = -fl qa smh9 q % ■ 
The displaced, orthonormal Fock basis constructed for each j a will be denoted 



la. 



(48b) 

and reads explicitly 
(49) 



As in the pure state framework, the set {a ■ \ . . .n q ^ ai rh q ^ a . . .)} provides an orthogonal operator basis for the to- 
tal lattice-phonon system, and the analogy can be carried further, step by step. The self-consistency conditions 
for the ansatz (EG) can be derived now from the expansion of the associated equation of motion (32) in the basis 



•)}• 



But, a word is in order first regarding the gauge generator K. One should note that substitution of the reciprocals 
of transformations ( fTi] ) in the bare phonon Hamiltonian H p h = ^huj q h q h q produces nondiagonal products of B's 

q 

and B's , which may prove cumbersome. But if it is observed, as in TFD jl2), that 

h\h q - h\h q = B\B q - BjB g , 



the gauge generator can be suitably rewritten 

K^K-^^bJb, , 



(50) 



which brings the free phonon Hamiltonian for Eq. ( |32| ) to the invariant form 

H ph - Hph = ^2 huj q ( b 9 b <? _ bjb 
q 

The equation of motion ( |32] ) acquires thus the particular form 

Hiat + w + (n ph - n ph 



in—— 

dt 



K 



7 • 



Subsequent substitution of the product ansatz (|40|) produces 



da , 
Hi 



7a 



ih 
~2 



(51) 



(52) 



a: , q 



+i 



Ti*^2/3 qa cosh^ g d • | lg ta ) - ih^^Pqa sinh# 9 d • | l q . a ) = 
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a | & * 'J a i 

a a a,Q 

+ hCJqPqa COsh 6" g a • | l gjQ ) + ^qPqa sinh ' I l<?,a) > (53) 

which is to be contracted, succesively, with each of the basis operators {a ■ \ ... riq tCn m q a . . .)}. 

The results of the preceding section can be extended straightforwardly, mostly by obvious substitutions. For 
instance, the contraction with operators {&' • | . . .n q . a > ,rh q , a > . . .)} carrying more than one excited 'quantum', i.e., 

^2 (n q + rh q ) > 1, shows that the interaction-gauge term I W — K J must be such that 
q ^ ' 

(...n g ,a,m,,a...|W-K|7 a ) = 0. (54) 

Hence the allowed phonon dependencies are zero order, first order or bilinear of the type B^B g ', B^B g ' and the 

tilde conjugates. However, it can be checked that, due to the specific character of W and K [as functionals of 
right-acting and left-acting operators], the bilinear terms can only arise accompanied by prohibited terms of the 

form fcosh^g sinh0 g /BjBj,^ and ^cosh# g sinh^q'BgBq'^ . For this reason, the allowed expressions for W and K at 

finite-temperatures [sinh q ^ 0] reduce to the very simple forms [in terms of zero temperature operators] 

W = £(w g bt+wtb 9 ) (55) 

and, respectively, 

K = K lat + J2 (v,b+ + v q bq) , (56) 



where the zero-order term in W is assumed subsumed in Hi a t, and K; at , w q and w q are lattice operators, functionals 
of c n and <y n only. At zero temperature, bilinear coupling terms can still coexist and the outcome parallels the results 
of Sec. II. Let us consider here the finite-temperature case. 
From the same contraction procedure it also follows that 

(a' | (1,, Q |W - K| %) | a) = {&' | (l q . a |W - K| %) \&) = Q, Vg, a £ a' (57) 

and 

^ = (&' \H lat \ &) + (&' \(%\W- K\%)\ a), Va^a', (a | 7) + . (58) 

Using expressions fl55|) and (|5^), one obtains from Eqs.(j57|) 

(a' |w g | d)cosh6' 9 — (a' |vt| a)sinh9 q = 0, Vq, a ^= a' , (a | 7) ^ , (59a) 

(&' a) sinh 0, - (a' \v q \ a)cosh<9 9 = 0, V<?, a f a', (a | 7) ^ , (59b) 

and similarly, 

( 7 a |W - K| %) = £ (/%*w g + /SgawJ) - K lat , (60) 
1 

which, inserted in Eq. (|58j), leads to 
da 



%h I a' 



dt 



(a' |H ex - K lat + + (3 qa ™\) I a) , V a £ a', (a | 7) ^ . (61) 
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As before, since a' spans a complete lattice basis, Eqs.(59) and ( pl| ) indicate that contributing d's, (d | 7) 7^ 0, satisfy 
operatorial equations of the form 



cosh#g Wq • d — sinh^q vj • d = rj q . a a 



(62a) 



sinh# g • d — cosh 6^ • d = 7y 9:Q 



and, respectively, 



in—— — 

dt 



(62b) 



(63) 



where ^ g . Q , fj q a and f2 Q are scalars. The explicit expressions for the latter can be identified from the contractions of 
Eq. ( [53| ) with the basis operators a' ■ \ lq\ a ')> d' 1 1 lg',a')) an d, correspondingly, a' • x f a i. A straightforward calculation 
yields 



Vq 



a = (& I (l 9 , a |W - K| 7 a ) | a) = U^jSga - TlLO q j3 q( }j cosh9 q , 



(64a) 



and 



■q q . a = (d I (lg, Q |W - K| %) I a) = (ih/3 qa - Tiuj q (3 qc }j sinh0 9 , 



(64b) 



i7 Q = ih d 



^ ) -(n H,„, n! + (n |,-,, \\ K -„S|n): 



£ fiw, |/3 ?Q | 2 - ^£ (/3 9Q /3* Q - p q J* qa 



(65) 



After substituting expressions (|64j), Eqs.(|62|) can be rearranged in such a manner as to separate the contributions 
in w q from those in v q . To this end, multiply Eq. (|6^a) on the right by (d^ cosh# g ), and the adjoint of Eq. (p2|b) on 
the left by (— d + sinh# 9 ), and add to obtain 

cosh 2 q w q ■ (dd^) — sinh 2 9 q (dd^) ■ w q — (ififl qa — huj q (3 qo ^j (dd^) . 

Similarly, multiply Eq. (|62|a) on the left by {a) sinh(9 g ), and the adjoint of Eq. (|6^b) on the right by (— d^ cosh6* g ), 
and add again to get 



cosh 2 9 q v q ■ (d^d) — sinh 2 9 q (d^d) • v q = 



(66) 



If we recall that d^d = Uaaa'U^ for some unitary transformation U a , and introduce the diagonal decompositions 
dd f = \Xk,o) v K, a (Xk,o.\ and « f « = J2 \x.K,a) v K ,a (x.K,a\, with \x,K,a) = U a \x K ,a), from Eq. (|flb) it follows that 



(cosh 2 d q v K > l0t - sinh 2 9 q v K . a ) (x K , a \v q \ Xk',<x) = , 



hence, 



(Xk,q \Vg\ XK',a) = {Xk,c \ v \\ Xk' ,a) = 

whenever v K>a ^ and/or v K i ;CX ^ 0. But then 



d • v q — a ■ v q = 



must hold, and Eqs.(|62|) are so reduced to the simpler form 

d = (iTif3 qa - fiUOqPqa ) & , 



(67) 
(68) 
(69) 
(70a) 
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w\ ■ a = 



(i%l3 qa - nw q /3 qa ) a . 



(70b) 



Equations (|69| ) show that all gauge couplings v q compatible with the ansatz ( (40[) for 7 give null contribution to the 
problem and cannot be employed for an eventual alleviation of computational complexity. Finally use of the quasi 
diagonal representation a = \Xk,o) \J v n,a (Xk.oI [corresponding to the diagonal representations for dd^ and d^d 

given above] in Eqs . (|70|) , leads to 



0, 



(71a) 



= 



(71b) 



I \Xn,a) ~ (ifrfiqa - hujgPqa^ \x K ,a) 

and reveals that the states \xn,a) contributing to a \v K ^ a 7^ 0] must necessarily satisfy the eigenvalue equations 

W g \Xk,o) = (iTlfiqa ~ huJqPqa^ \Xk,o) , (72a) 



v\ |X«,a) = (ift$qa ~ huj q l3 qo }j \Xn,a) 



Finally, in view of Eqs.([7G) and (pq), the effective evolution Eq. (p3|) for d reduces to the unperturbed form 



da 
in—- = 
dt 



H/ at — K; at — fl a 



(72b) 



(73) 



It is immediate that Eq. ( |73[ ) preserves the orthogonality of distinct d's, i.e., d(a'\a) /dt = 0, as required by the 
ansatz for 7. 

While the self-consistency constraints for the pure state case assumed the form of modified conditions for 
decoherence-free propagation, Eqs.(^) and (|7^) show that in the finite-temperature case, self-consistency corresponds 
exactly to the conditions for decoherence-free decoupling. The reason for this is seen in the effect of the interaction 

term — Jij on any D 2 product, i.e., 



(w-k) 



a ■ 7 Q 



K iQt d 



(74) 



which means that the interaction becomes decoupled, regardless of the exact nature of the state for the phonon sub- 
system. Here, denotes the eigenvalue of w q corresponding to a [^ q , a — Vq,a/ cosh. 9 q ]. 

We may conclude now that a finite-temperature state operator of the form ( f40| ) is an exact solution of the evolution 
equation fl52| ) if and only if 

a) the lattice-phonon interaction W is linear in the phonon degrees of freedom [Eq. (p5l)1 ; 

b) the lattice operators d are state operators on decoherence-free subspaces of the lattice subsystem, i.e., on [degen- 
erate] common eigenstates of w q and w q [Eqs.((70]) and ([72])], which are, simultaneously, solutions of the evolution 
equation (|73|); 

c) the phonon displacement parameters (3 qa evolve according to 



iTlj3q a TliO q [3 qa — f^q,a 



(75) 



In addition, 

d) compatible gauge generators act on the lattice states only, Kd = K; at d for all d, and the lattice energy shift fl a 
in Eq. (f73|) amounts to [see Eq. (EE 



n a (t) 



^^2(Vq a f3 



(76) 



Obviously, this result is strongly reminiscent of the pure state case examined in Sec. II, but involves a simpler form 
of W and the additional constraint (|72jb) on the lattice states, which calls for exact decoherence free propagation in 
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the lattice subsystem. It should be noted, however, that this constraint vanishes in the zero temperature limit, when 
sinh0 9 = 0, and the pure state case of Sec. II, for linear phonon coupling, can be recovered identically if needed. For 
the particular situation when the initial phonon displacements are null, j3 qa = 0, and the initial state is a product 



between a lattice distribution and a phonon thermal state, p (0) 
theorem, referred to in the introductory section: 



Plat (0) p P h,T, one obtains the following interesting 



A system initially in a product state p (0) = pi a t (0) Pph.T evolves into a Gaussian, generalized D\ state given by the 
ansatz if and only if the interaction is linear in the phonon coordinates and the initial lattice state is a distribution 
on a direct sum of orthogonal decoherence-free subspaces, provided any exist. 



When the eigenvalues p q . a are time independent, Eq. (ji 
harmonic oscillations around displaced equilibrium positions, 



5l) shows that the displacements /3 ga perform simple 



qa (t) 



P qa (0) 



TlUJ n 



(77) 



In the peculiar case when the initial displacement brings each mode over the displaced equilibrium position, the 
reduced phonon density-matrix becomes time independent [see Eq. (^4|), and take into account that (d|d) = const., 
according to Eq. fl73|)]. But the stationary phonon state is not exactly thermal, unless the overall state is a simple D2 
product, and the lattice evolution is confined to a single DFS. Furthermore, since all generalized amplitudes (-y a \ j a i) 
are constant in time, the lattice reduced density-matrix pi at = (7a I 7a') d'd^ evolves in an unperturbed fashion, 

a, a' 

up to the [constant] energy shifts Q a , i.e., 



Plat 



-jr(H lat 



a (0) d f (0) exp 



h {Hlat 



As for pure states, the evolution of self-consistent, density-matrix Di states ( p0| ) is disentangled into separable 
evolutions of the component D2 states d • j a , hence no self-consistent model can generate equations of motion 
coupling distinct D2 contributions. Again, no nonlinearity survives in the equation of motion for the lattice states, 
and standard Davydov distributions, for which d = <Pra c n |0)/ a t (Xa\, for some orthonormal \x a ) [{Xa I Xa>) — $aa']> 
can generate, when self-consistent, only stationary lattice configurations [(d|d) = \ f n \ 2 — const.]. A similar statement 
applies to multiquanta, on-site D\ states and their superpositions. Also, in complete analogy to the pure state case, 
the relation to the D 2 Davydov soliton equations becomes apparent when the equations of motion for d and f3 qa are 
cast in the form 



da 
in—— = 
dt 



Hi Qt - Ki a t + (7a |W| la ) + n a 



(78a) 



ih/3q a - TlLO q (3 qa = (& \w q \ d) 



(78b) 



If d is given a linear parametrization, one can recover Davydov-like equations, but the self-consistency constraints 
([72|) for d reduce the interaction term in Eq. (|78|a) to a scalar contribution. Even this scalar cannot have a nonlinear 
dependence on d, since the right-hand side of Eq. (f78|b) becomes independent of d. 

At last, the self-consistency conditions (^21) can be stated alternatively as the initial state constraint 



A Q d (0) = 



(79) 



with 



00 

A a = J dt£ [(wt jQ (t) - Pq a (t)) ( Wq>a (t) - /i q , Q (t)) - 



+ (w g>Q (t) - p q , a (t)) (wj (t) - fj* ja (t))] 



(80) 



and 
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w g , a (t)=U+ (t)w, (i)U Q (i) , 



(81) 



where U Q (t) is the unitary propagator corresponding to the equation of motion ( |73| ) . 

Recall that the above results apply to the finite-temperature version of the ansatz (|i(]) . It is interesting to keep in 
mind that the zero temperature limit also provides a density-matrix generalization of the pure D\ states, for which the 
acceptable lattice-phonon couplings resume expression (pd|), bilinear in phonon coordinates. This allows, conceivably, 
an extended class of nontrivial models and the problem can be approached in the same manner. Nevertheless, only 
models with interactions linear in the phonon coordinates admit self-consistent, finite temperature extensions, and 
thus bear realistic physical significance. One may note that both the original Davydov model Jig] and all its studied 
versions |Q belong to this class. 



V. DISCUSSION AND EXAMPLES 



For the reasons outlined above, we limit the discussion to models with linear phonon coupling, and extend the 
self-consistency conditions for pure lattice states to include the constraint 



wt la) 



^qPqa 



a 



(82) 



and place the state \a) in a DFS of the lattice. In other words, \a) can be a self-consistent D2 factor if and only if it 
lies in a lattice DFS. The properties of the decoherence-free subspaces have been studied in detail in connection with 
quantum computation theory, and the interested reader is referred to the available literature ^fft- Suffice it to say 
that their existence is determined by the symmetry properties of both the unperturbed Hamiltonian and the system- 
bath interaction, and thus, a proper understanding involves a Lie-algebraic/group-theoretic framework. However, 
for an expedient assessment of various Davydov-like models, the brief pedestrian characterization given below proves 
satisfactory. 

Let us start by noting that instead of Eq. (p^), it is sufficient to require that |a), which according to Eqs.(|22|) and 

(p^), satisfies w q \a) — (ih$ qa — TiWqPqaj \ a )> must belong to an eigenspace S of w q that is also left invariant by w q . 

Indeed, since w q reduces to the identity on S, up to a scalar factor, if w q leaves this subspace invariant, it also must 
be proportional to the identity, by the complex conjugate factor, and the adjoint eigenvalue equation ( |82] ) is recovered 
necessarily. Since \a) belongs to a DFS, the equation of motion for |a) reduces to the unperturbed form [see Eqs.(^6[), 
@, and ©] 



ih Id) 



lat 



(83) 



Hence, in all cases that support a smooth extension to finite-temperature conditions, the self-consistent lattice states 
are effectively decoupled from the phonon dynamics, up to a phonon modulated phase factor. 

In fact, the remarkable phenomenon highlighted by decoherence-free propagation is that coupling to a [nonequi- 
librium] thermal bath need not necessarily result in decoherence to a mixt, statistical state. In particular, any 
self-consistent D2 product with a pure lattice factor \a) displays a finite temperature density-matrix 



p (t) = exp 



-Hi at t 



\a(0))p P h ({P ga },T) (a (0)| exp 



-Hi at t 



with 



P P h ({P qa },T) = l[-^ 



exp 



TlUJ q 



H - P* qa (t)] [b q - f3 ga (t)] 



It obviously maintains the pure character of the lattice state throughout the evolution, provided the time-dependence 
of the displacement parameters complies with the self-consistency conditions. For a simple concrete example, let the 
lattice system be coupled to the phonon bath in a site-homogeneous manner, such that, e.g., w q — \q X) c n c n = Xa-^laU 



and let the unperturbed lattice Hamiltonian conserve the number of lattice excitations, i.e., 



W q ,Niat 



w\,Ni at 



Hlat , Ni a t 



= 0. Then any lattice state with a well-defined number of excited quanta qualifies as a D2 factor state, 



and for Ni a t |a) = v |a), the associated equations of motion for the displacement parameters read iTi(3 q — Tiuj q f3 q = v\ qi 
with the trivial solution (3 qa (t) = \(3 qa (0) + {vxq/fi^q)] exp(—iuj q t) — (vxq/huq)- Hence, although the lattice interacts 
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with a [thermal] bath, its state is propagated according to the unperturbed dynamics, and if the original state is a pure 
state, its condition is preserved in time. On the other hand, the bath evolves into a superpositon of Gaussian states, 
even when initially in a thermal equilibrium state [f3 qa (0) — 0]. Only if at the outset the bath modes are displaced 
directly over the displaced equilibrium positions (—isx q /huj q ) , does the bath remain in a stationary state. The effect 
is absolutely robust under variations of the interaction strength, i.e., stronger coupling cannot induce decoherence. 

Returning now to the formal characterization of decoherence-free propagation, observe that if the lattice couplings 
w q and w q are time independent, \a) must dwell in a time-invariant subspace [corresponding to time independent 
eigenvalues fi qa ] and, since H\_ a t alone drives the evolution of |a), this subspace must also be an invariant subspace of 
Hi a t [the latter can be time dependent]. To see this, let us assume that the eigenvalue p, qa corresponding to \ct{t)) 
changes parametrically in time, such that the time derivative of w q \a) = [i qa |a) reads w q \a) = fi qa \a) + (dfi qa /dt) \a). 
If \ct±) denotes the component of |d) which is orthogonal to all eigenvectors of w q [wj] for the eigenvalue fi qa [Hq a ], 
then Wq \a±) = /i qa \ct±) + (d^i qa /dt) \a). But the latter equation implies that (a\ w q \a±) = (dn qa /dt) (a \ a) = 
[since (a\w q — [i qa (a\], and thus {d(i qa /dt) = 0. In fact, the requirement that self-consistent states |a) must lie 
in an invariant subspace of the lattice Hamiltonian extends to all cases where the couplings w q have time-invariant 
eigenspaces, but not necessarily time independent eigenvalues [e.g., when w q — J^Xqn {t) c n c n\- Applying the same 

n 

reasoning, it follows that in such situations the time-dependence of the eigenvalues corresponding to a self-consistent 
| a) is limited to the intrinsic time-dependence implied by the form of coupling, i.e., (dfi qa /dt) — (dfi qa /dt). In the 
most restrictive situation, when \a) is a nondegenerate eigenvector for at least one w q [w q ], then it can be solution of 
Eq. (|83|) if and only if it is also a time independent eigenstate of Hi at . 

The necessary and sufficient condition that a common eigenspace S of all w q and w q be a decoherence free subspace, 
i.e., that it remain a common eigenspace of the w q 's and wj's under the dynamics driven by the unperturbed lattice 
Hamiltonian, is that every \a) in S must also satisfy 

('[„,.,, „, ] + ^) W0) = ^ Wi)) , (84a) 



dfJ^o 
dt 



\a(t)) , 



(84b) 



for every q, at all times. The proof is trivial. Let |a(£)) evolve according to Eq. ( |83| ) and satisfy w q (t) \cx(t)) = 
Mg,a (*) \ a (*)) a * some instant t. If w q (t + At) \a (t + At)) = \i q . a (t + At) \a (t + At)) also holds, then one has, to 
first order in At, 



w q (t) \a (t)) + At 



w q (t)\a(t)) + ^\a(t)) 



(t)\a(t))+At 



Vq,a (t) \a{t)) 



which in view of Eq. ( J83| ) yields Eq. (p4a). Conversely, if \a (t)) also satisfies Eq. 
slight rearrangement of terms gives 

' l w q (t)\a(t))} = f t [fx q , a (t)\a(t))} , 



) at the same instant t, then a 



dt 



and Wq (t + At) la (t + At)) = [i q a (t + At) \a (t + At)) follows necessarily. A similar reasoning applied to w q verifies 
condition Eq. (|84[b). If S is time- invariant, then the restriction of w q to S changes in time only by a time-dependent 
scalar multiplication, w q (t) \s — X (t) w q (0) \s, and 

^l«(t)> = ^pi«(t)> , 
such that we obtain the following restricted criterion: 

A time-invariant common eigenspace S of all w q is a DFS, i.e., is left invariant by both Hi at and every w q , if and 
only if it is contained in the common kernel of all commutators [w q , w q ~\ and [Hi at , w q \. In that case S is also in the 
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kernel of every [H lat ,w\] [i.e., w q ,w q , |a) = [H lat ,w q ] |a) = [Hi a t,w q ] |a) = for every q, and every \a) in S]. 
The above can be given a proof independent of Eqs.(|84|). Indeed, assuming that S is time independent, if w q and 



a 



Hi a t leave S invariant, for any \a) in S it is true that w q w q \a) = [i qa w q \a) = w q w q \a) and similarly, Hi at w, 
HqaHiat \ct) = WqHi a t \ot) • It is also true that, for all |a) in S, (f3 \w q \ a) = whenever ((3 \ a) — 0, and (a' \w q \ a) = 
Mg Q («' I a ) f° r au \ a ) and K) m S. Hence |a) = fj,* a \a) for any \a) in S. But then Hi at w q \a) = fi* a Hi at \a) = 
w\Hiat | a)) because Hz Q t |a) is necessarily in S. Conversely, if [u> g ,u>J] \a) = for all \a) in S, then vA \a) is in S, 
(a' \w q \ a) = fi* a (a' \ a) for all \a') in S and w q \a) = fi* a \a). And if [Hi at ,w q ] \a) = for any \a) in S, then H\ at \a) 
is an eigenvector of w q for the eigenvalue fJ, qa , and must be in S. It is also immediate that [H/ a t,wJ] |a) = holds 
too. 

Virtually all studied versions of the Davydov model display time independent interactions, and therefore fall under 
the incidence of this prescription for self-consistency. For example, the general Frohlich Hamiltonian 

H = ^ J-mnclnCn + ^ TlU q b\b q + ^ TlW q (Xqn,b\ + X qn b q) c l c n (85) 
mn 

employs the couplings 



q qn 



1 </ 



= hbJ q Vx g „4 C « ) V <?> ( 86 ) 



which obviously satisfy [tu g , toj] = and have as eigenvectors the unperturbed on-site Fock states, for the discrete 
eigenvalues fi qn = TiuJqXqn- If all \ qn are nondegenerate for at least some q, as happens in the proper Davydov model, 
where Xqn = Xq etnaq j the lattice Hamiltonians admitting self-consistent D2 states under the interactions ( J86| ) can 
only be of the diagonal form Hi at — '^2 i e n c\ l c n . This is, of course, the well-known result that the D\ states are exact 



under the phonon-lattice interaction ( J86[ ) just in the limit of vanishing hopping between monomers. Only when the 
Xqn collapse to degenerate values, e.g., in the low-frequency limit, when qa — > 0, can nontrivial D2 states appear for 
the Hamiltonian (|85|). But in order to recover a nontrivial dynamics in the low- frequency limit, the phonon modes 
must belong to an optical band, whereas the Davydov model uses specifically acoustic modes. 
The situation does not improve when hopping terms are added to expression (|86|), to obtain 

W q = flUJ q ^2 [Xqllclcn + Xq^ ( c l+l C n + c Un+l 



\xq 1] CnCn + X q 2) (4+l C n + 4c«+l) 



a maq 



Vq, 



(87) 



as done by Takeno (T(| , Todorovic et al. , Bartnik et al. or Pang . To the contrary, in this case w q and w q 
no longer commute, and do not share eigenvectors, which compromises self-consistency from the outset. For instance, 
with the translation invariant ansatz in the second line above, one obtains 



[w q ,w, 



.(2)' 



'n-l C «+l 



■ 1^71—1 



(88) 



with 



21m 



1) ( (2)Y 
[Xq ') 



and 



(2) _ 



-2 



(2) 
Xq 



sin (qa). But although [t 



..tl 



is diagonalized by 



the simple Fourier transformation c„ 



Wq,W q \ 

-jjj ^2 CQe l ® na and can have null eigenvalues, none of its kernel states is 



an eigenstate of w q or w' q . As for the proper Davydov model, self-consistent D2 states can only appear in the 
low-frequency limit, provided the phonons belong to optical modes. 

In view of the above, a search for improved wave functions has no chance of recovering a self-consistent model 
as long as the structure of the Hamiltonian is not adjusted to support decoherence-free D 2 states. The partially 
dressed state, introduced by Brown and Ivic ]20[ ], which belongs yet to the standard D\ class, can yield a better 
approximation for the exact dynamics of an initial D\ configuration under the Davydov Hamiltonian, but does not 



make any progress toward a self-consistent model. The same holds for the coherent state used by Wang et 



their vibron soliton model, and for Pang's quasicoherent wave function |<3>) = i 



1 + 



4 



21[ in 
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|L9|. In these cases it suffices, in fact, to note that the interactions ( |8q ) and ( |87| ) commute with the lattice number 
operator, N~i at = ^ c n c n, hence, any eigenvectors are also eigenvectors of Ni a t- But a coherent or quasicoherent state 

n 

is not an eigenvector of Niat, and cannot generate a self-consistent D2 state even in the low- frequency limit. Yet the 
underlying idea, that the lattice state \a) could be a coherent state [or an approximation thereof], does provide a 
good starting point for the construction of a nontrivial self-consistent example. 

For instance, consider the widely used version of the [time independent] Frohlich Hamiltonian for a one-dimensional 
lattice of bosonic oscillators with nearest-neighbor hopping, 



H = ^2 ec « c » ~ J Y ( c i+i c » + 4 c «+i) + huj q b l b q + huj q (x q nb\ + X* qn b q) 4i°n , 



Ql) 



to which we add an external pumping term of the form 

bq) (Cn (*) cj, + C (*) Cn) 



(89) 



(90) 



nq 



In such a case, the corresponding couplings w g are readily identified as 



W q = huig 



£ **» ( C « + Ct (*)) (Cn + Cn (*))-£ X 9 n |Cn (*) f 



(91) 



satisfy [u> 9 , = 0, and have obvious time-dependent eigenstates. Let us consider the ground state 



e « e "(') exp 



-X)(C»(i)4-CWc„) 



10) 



Since |a) must be simultaneously a solution of the unperturbed problem 

d\a) 



ih- 



dt 



Y £C l C n ( C n+l C ™ + c nC«+l J - 

mn n 

it is necessarily a product of coherent states for the lattice normal modes , with a phase factor 

e a (t) = - y d T n a ( T ) . 

Hence the pump-induced displacements must amount to Cn (i) = X) Cg (0) e mqa e~ lUqt 1 and satisfy 



(92) 



(93) 



(94) 



ih —JT = £ Cn ~ J (Cra-1 + Cn+l) 



(95) 



where the lattice frequencies read ?if2,j = e — 2 J cos (qa). The corresponding equations of motion for the phonon 
displacements are now 



ih 



dt 



fiUJn 



(96) 



and apply at both zero and finite-temperature. This is a simple example showing that external pumping can be 
instrumental in maintaining the lattice in a pure coherent state at arbitrary temperature. Moreover, it also verifies 
that a Frohlich-like pumping effect can be produced, albeit on the bath modes, due to the nonlinearity in Eq. J96|). 
Indeed, let us assume, as in the Davydov model, that the bath phonons belong to an accoustic branch of N distinct 
modes, such that \qn = Xq e ~ lqna an d u q = oj- q . Taking advantage of the degeneracy of the lattice spectrum, let the 
external pump act at the lattice frequency Qq, such that Cn = (£ + e l ® na + Q_ e ~ l Q na ^ e - in Q t , But then the driving 
term for the phonon displacements in Eq. (Bq) amounts to 



Xq J2 e ~ iqna IC" Wl 2 = N M [flC+l 2 + IC-I 2 ) S q .o + C+C-^,2Q + C+C-S q .-2Q 
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which means that only the fundamental bath mode and the modes of frequency uj 2 q experience time independent 
external driving. Therefore, if the phonon bath is initially in thermal equilibrium, all modes that are not affected by 
the external pump remain in thermal equilibrium, while the three modes q = 0, ±2Q are each driven in a Gaussian 
state, eventually of macroscopic displacement [see the presence of the N factor above] . In another limit case, when 
the initial displacements of the bath modes coincide with the displaced equilibrium positions set by the amplitude 
of the pump, the bath remains in a stationary, nonequilibrium state, while the lattice oscillates harmonically at the 
pump frequency. Unless Q — or an accidental degeneracy intervenes, the pumping frequency for the coherent lattice 
state is distinct from that of the coherently driven modes of the bath, and the process obviously qualifies as a type 
of Frohlich effect. But contrary to the usual picture, it is the bath that is driven into a nontrivial state, while the 
primary, lattice subsystem experiences no nonlinearities, and is maintained in a pure state. 

For closure, let us illustrate more fundamental features of the coherent-product ansatz [Eq. (|40|)] in another simple 
example with exactly tractable decoherence-free propagation. Consider a symmetrical dimer with site-independent 
Frohlich interaction and phonon-assisted tunneling, described by 



H = 



(c\c± + c\c 2 
The simple transformation 



(cl 



c 2 + c\c\ 



c 3 = ±= (cx + (-ir ^ 



.1 = 1, 2 



brings the Hamiltonian (|97j) to the evidently solvable form 

H = £ic\c x + £ 2 c\c 2 + c\c! 2J (Xq,l b q + X*q.l b q) + 4 5 2 2j U<Z,2&J + Xq,2 b q) + X! ^9^9 ' 



(97) 



(98) 



(99) 



where the renormalized energies and coupling constants read £\ = s — J, £ 2 = £ + J, Xq.i — Xq ~ \ and, respectively, 

Xq,2 = Xq + \- Since H, c\ci = H, c 2 c 2 = 0, the interaction with the bath does not mediate an energy exchange 

between the monomers, but affects only their coherent correlations through a process of dephasing. In its fermion/two- 
level realization, this model is easily recognized as equivalent to the exactly solvable Jaynes-Cummings model of 
quantum optics [^2| and to the Caldeira-Leggett model of quantum dissipation j23j. Various versions have also 
enjoyed attention lately in studies of decoherence in quantum registers p^Ps| - Confining the discussion to single- 
quantum lattice states, so as to cover both the fermionic and the bosonic case, we apply the ansatz ( |40| ) and construct 
the self-consistent, statistical superposition 



7 = ¥>l(f) 4\0)lat(<t>l\ 71 + 4\0)latifa\ 



72 



(100) 



with ifj (t) = <fj (0) exp 
lattice state vectors, and 



scalar amplitudes, 



7j = ]I ex P [P q] (t)b\-P* q] {t)b q 



normalized, but not necessarily orthogonal, 



: exp 



2k B T^ bq 



S 2 ^ 2 -. The phonon driven energy shifts flj are given by Eq. 



where (3 q j (t) 

appropriate substitutions. Then, the total density-matrix for the dimer-bath system reads 



p = |1> (1| 7i7i f + \<P2\ A |2) (2| 7 27 2 T + Vl<P*2 (0i I 02) |1) (2| 7 i7 2 T + V*i^2 (02 | 0i) |2) (1| 727! , 



(101) 
|) with the 
(102) 



where we have denoted, for simplicity, |1) — c\ |0), . and \2) — c 2 |0), and the reduced density-matrix for the dimer 
is obtained accordingly as 



Tr ph p =|^| 2 |I)(I| + |^ 2 | 2 |2) pl + tptrt 



02) (72 I 7i) |1) (2| + ^1^2 (02 I 0l) (71 I 72) |2) (1| • (103) 



Since \ipj (t)\ = \tpj (0)|, the occupation probabilities of the states |1) and \2) are invariant, and only the correlation 
(/>dimcr)i2 = (^dimer)lj = <Pi<p2 (0i I 02) {l2 \ 7i ) varies in time, as expected in a typical dephasing process. Of course, 
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the same is not true for the original, localized states |1) = c{ |0) ia( = (|2) — |1)) and |2) = |0) tet = 4= (|1) + |2)). 
In terms of the latter, the dimer reduced density-matrix reads 

Pdimcr = Q + R-e (pdimer)l2^ I 1 ) (1| + Q ~ R« (Pdimer)i2^ I 2 ) < 2 I + 



(\h\ 2 ~ 1^1 1") +ilm(p d imer) 15 ) |1) (2| + h.C. 



(104) 



and it becomes apparent that the correlation between the physical dimers disappears if and only if \4>i\ 2 = \4>2\ 2 = \ 
and Im (pdimcr)i2 = 0. It is also evident that all nontrivial aspects of the dynamics are carried by the matrix element 
(/5dimer)i2) an d particularly by the bath factor 



(72 I 71) = (7r |exp 



J2(P q 2 (t)bl-P* q2 (t)h q ) 



■ exp 



J2(P q i (*)bt-/3* x (t)b q ) 



It) 



(105) 



Expression ( |105| ) can be evaluated straightforwardly by using the thermal operators (J47|), and amounts to (72 | 71) 
exp[-i$(i)] exp [-r (*)] , where $(t) = i £ (/3,i/3* 2 - and 



r(t) 



1 v 2 *^/ 



(106) 



Under the traditional assumption of initial thermal equilibrium of the bath, with null original displacements, the 
initial density-matrix corresponding to the ansatz (102) factorizes in the usual manner, as p (0) = pi at (0) <g> p p h,T, 
where 



Plat (0) = |V1 (0)| 2 |1}(T| + |^ 



|2) (2| + <^ (0) ^ (0) <0r I 2 ) |1) <2| + <p\ (0) ^2 (0) (02 I 0i) |2) (1| , (107) 

1 — l) , the relaxation 



and p p h,T = ItIt- Since at later times the bath displacements read f3 q j (t) = {Xqj/^g) ( e 
function T becomes 



r° (t) =2^|A g | 2 coth 



TlUJ q \ (1 — COS UJqt) 



2k B T 



(huj q 



(108) 



The physics behind expression (108) has been discussed at length in Rcf. |25[. Depending on the spectral density 
of the bath modes, the relaxation exponent may or may not have a finite limit as time passes to infinity. If it does 
not reach a plateau, the correlation (pdimer)j2 eventually falls to zero and the dimer density-matrix is invariably 
driven toward a steady distribution diagonal on the states |1) and |2), hence, undergoing total dephasing through 
an environment-induced superselection process pj] ]. The noteworthy point in the above derivation is that the dimer 
density-matrix (107), parametrized by three linearly independent, real parameters [\(j>i\ 2 + |<fe| 2 = 1], may spa n the 
entire set of statistical dimer states corresponding to the subspace {|1) , |2)} [or {|1) , |2)}]. In other words, Eq. (102) 
gives explicitly the general solution of the Liouville-von Neumann problem for the Hamiltonian (|97]), with an initial 
condition of the type p (0) = pi a t (0) ® p P h,T- 

From the point of view of coherent behavior, a special mention is reserved once more for the singular situation 
when the bath displacements remain stationary in their displaced equilibrium positions, i.e., (3 q j (t) = — (x a j / fru> ) . 
Then the bath is left in a steady superposition of thermal Gaussian states, and the corresponding bath factor ( |l06| ) 
becomes constant in time. If the latter does not vanish, the matrix element (/3dimcr)i2 oscillates as (pi (t) (p% (t) = 
(pi (0) ifi2 (0) exp [— iAfii] , at the frequency 



ill 



1 



hu>„ 



(109) 



and the dimer is kept in a coherently oscillating [or two-level rotating] state, somewhat reminiscent of a soliton state. 
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VI. CONCLUSION 



We have analyzed the self-consistency conditions for a fairly wide class of generalized D\ states, based on the 
Davydov ansatz for soliton propagation in molecular chains. Our extended ansatz is given in a density-matrix [state 
operator] form [Eq. (f4(i|)], which covers, in proper limit cases, both the pure state and the finite-temperature standard 
Davydov states. In general, it describes strongly entangled, statistical superpositions of lattice and bath states, with 
the property that the bath is maintained in a statistical superposition of Gaussian thermal states. We find that the 
exact propagation of such states amounts to an effective decoupling of the lattice subsystem from its boson bath, i.e., 
to decoherence-free propagation. Given the specific form of the bath state, the lattice-bath interactions compatible 
with such a phenomenon are limited to linear or bilinear forms in the bath coordinates, but only interactions linear in 
the bath degrees of freedom allow self-consistent propagation at finite-temperatures. Unfortunately, in all cases, the 
equations of motion for the ansatz parameters differ from Davydov's soliton equations by additional constraints, which 
apparently abolish the characteristic nonlinearities of the latter. Just as Brown's validity theorems j|, the result is 
independent of any internal characteristics of the system, including symmetries, range and strength of interactions, 
or type of phonon branches. This does not affect the suitability of the ansatz as a variational trial ansatz, and the 
Davydov model may still be a good approximation in different conditions. 

As an interesting byproduct, we are left with a set of nontrivial, exact density-matrix solutions for systems sup- 
porting decoherence-free propagation. Although exact results for such systems have been already reported [see, 
e.g., |ic|]25j] ], the present solutions are given in a closed, explicit operatorial form, and are not limited to the usual 
equilibrium, separable, initial conditions. Consequently, it becomes possible to probe the dynamics of such systems 
under nonequilibrium, Gaussian states of the bath. We retrieved a Frohlich effect under decoherence-free conditions 
facilitated by time-dependent, external pumping, and found that certain nonequilibrium states of the bath can be 
instrumental in maintaining bath-entangled, decoherence-free states of the driven (lattice) subsystem. We believe that 
this ansatz may prove useful in illustrating a variety of nontrivial aspects of relaxation and decoherence in quantum 
systems. 



[1] D. W. Brown, K. Lindenberg, and B. West, Phys. Rev. A 33, 4104 (1986); B. S. Tosic, Lj. D. Maskovic, J. P. Sertrajcic 

and D. Lj. Mirjanic, Physica A 52, 217 (1988). 
[2] D. W. Brown, Phys.Rev. A 37, 5010 (1989); 

[3] For work before 1991, see the now classic review by A. C. Scott, Phys. Rep. 217, 1 (1992) and refs. therein; also "Davydov's 
Soliton Revisited" (NATO ASI Series B, vol.243), eds. P. L. Christiansen and A. C. Scott (1991); following is a short 
selection of more recent publications. J. W. Schweitzer, Phys. Rev. A 45, 8914 (1992); V. May and M. Schreiber, Phys. 
Rev. A 46, 6925 (1992); W. Forner, Theor. Chem. Acc. 101, 223 (1993); J. Phys.: Cond. Matter 5, 803 (1993); 5, 3897 

(1993) ; Phys. Rev. B 53, 6291 (1996); Int. J. Quant. Chem. 64 , 351 (1997); also Ref. @ below; A. R. Vasconcellos and 
R. Luzzi, Phys. Rev. E 48, 2246 (1993); L. Cruzeiro-Hansson, Phys. Rev. A 45, 4111 (1992); Phys. Rev. Lett. 73, 2927 

(1994) ; Phys. Lett. A 223, 383 (1996); 249, 465 (1998); V. M. Kenkre, S. Raghavan and L. Cruzeiro-Hansson, Phys. Rev. 
B 49, 9511 (1994); L. Cruzeiro-Hansson and S. Takeno, Phys. Rev. E 56, 894 (1997); Z. Ivic, Physica D 113, 218 (1998); 
A. M. Clogston, H. K. McDowell, P. Tsai and J. Hanssen, Phys. Rev. E 58, 6407 (1998); Xiao Yi, Phys. Lett. A 243, 1747 
(1998); J. Tekic, Z. Ivic, S. Zekovic and Z. Przulj, Phys. Rev. E 60, 821 (1999); Z. Ivic, Z. Przulj and D. Kostic, Phys. 
Rev. E 61, 6963 (2000). 

[4] A. C. Scott, Phys.Rev. A 26, 578 (1982); Phys. Scr. 29, 279 (1984); P. L. Christiansen, J. C. Eilbeck, V. Z. Enol'skii and 
Ju. B. Gaididei, Phys. Lett. A 166, 129 (1992); A. V. Zolotaryuk, K. H. Spatschek and O. Kluth, Phys.Rev. B 47, 7827 
(1993); A. V. Savin and A. V. Zolotaryuk, Physica D 68, 59 (1993); Z. Ivic, D. Kapor, M. Skrinjar and Z. Popovic, Phys. 
Rev. B 48, 3721 (1993); A. V. Zolotaryuk, K. H. Spatschek and A. V. Savin, Europhys. Lett. 31, 531 (1995); M. Daniel, 
and K. Deepamala, Physica A 221, 241 (1995); A. V. Zolotaryuk, K. H. Spatschek and A. V. Savin, Phys.Rev. B 54, 266 

(1996) ; Z. Ivic, D. Kostic, S. Przulj and D. Kapor, J.Phys.iCond.Mat. 9, 413 (1997); also refs. (r|- |l[ below. 

[5] Qing Zhang, V. Romero-Rochin and R. Silbey, Phys.Rev. A 38, 6409 (1988); M. J. Skrinjar, D. V. Kapor and S. D. 

Stojanovic, Phys. Rev A 38, 6402 (1988). 
[6] P. Zanardi an d M. Rassetti, Phy s. Rev. Lett. 79, 3306 (1997) [also as |quant-ph/9705044| ; Mod. Phys. Lett. B 25, 1085 

(1997) [also as |quant-ph/971004l| ]; D. A. Lidar, D. Bacon and K. B. Whaley, Phys. Rev. Lett. 82, 4556 (1999); D. Bacon, 
D. A. Lidar and K. B. Whaley, Phys. Rev. A. 60, 1944 (1999); S. De Filippo, Phys. Rev. A 62, 052307 (2000); D A. 



Lidar, D. Bacon, J. Kempe and K. B. W haley, Phys. Rev. A 63, 022306 (2001); 63, 022307 (2001); ^uant-ph/0004064 



Song and R. Szabo, ^uant-ph/0011021 



D. 



[7] L.-M. Duan and G.-C. Guo, Phys. Rev. Lett. 79, 1953 (1997); P. Zanardi, Phys. Rev. A 56, 4445 (1997); 63, 12301 (2001); 



21 



D. A. Lidar, I. L. Chuang and K. B. Whaley, Phys.Rev.Lett. 81, 2594 (1998); D. Bacon, J. Kempe, D. A. Lidar and K. 
B. Whaley, Phys. Rev. Lett. 85, 1758 (2000); A. Beige, D. Braun, B. Tregenna and P. L. Knight, Phys. Rev. Lett. 85, 
1762 (2000); P. Zanardi and F. Rossi, Phys. Rev. Lett. 81, 4752 (1998) : Physica B 272, 5 7 (1999); D. Kielpinski et al., 
quant.-ph/0102086; A. D. Greentree, S. G. Schirmer and A. I. Solomon, quant-ph/0103118; ; P. G. Kwiat, A. J. Berglund, 
J. B. Atlepeter and A. G. White, Science 290, 490 (2000). 

[8] M. Tegmark, and H. S. Shapiro, Phys. Rev. E 50, 2538 (199 4) [also as br-qc/9402026| . 

[9] M. Tegmark and L. Yeh, Physica A 202, 342 (1994) [also as |gr-qc/940100q |. 
[10] A. Venugopalan, Phys. Rev. A 61, 012102 (1999). 

[11] W. H. Zurek, Phys. Rev. D 24, 1516 (1981); Phys. Rev. D 26, 1862 (1982); J. P. Paz and W. H. Zurek, Environment- 
Induced Decoherence And The Transition From Quantum To Classical, vol. 72 of Les Houche s Summer School S eries, 
pg.533, edited by R. Kaiser, C. Westbrook and F. David (Springer- Verlag, Berlin, 2001) [also as quant-ph/0010011 . 

[12] H. Umezawa, H. Matsumoto and M. Tachiki, "Thermo Field Dynamics And Collective States" (Amsterdam, North Holland, 
1982); also Y. Takahashi and H. Umezawa, Collective Phenomena 2, 55 (1975); T. Arimitsu and H. Umezawa, Progr. Theor. 
Phys. 77, 32 (1987); 77, 53 (1987). 

[13] I. D. Lawrie, J. Phys. A: Math. Gen. 27, 1435 (1994). 

[14] W. Forner, J. Phys.: Cond. Mattter 4, 1915 (1992); 5, 823 (1993); 5, 3883 (1993). 

[15] A. S. Davydov and N. I. Kisluka, Phys. Stat. Sol.(b) 59, 465 (1973); A. S. Davydov, Zh. Eksp. Theor. Fiz. 78, 789 (1980) 

[Sov. Phys. JETP 51, 397 (1980)]. 
[16] S. Takeno, Progr. Theor. Phys. 71, 395 (1984); 73, 853 (1985); J. Phys. Soc. Japan 59, 3127 (1991). 
[17] D. Todorovic, Lj. Ristovski and B. S. Tosic, Phys. Stat. Sol. B 190, 251 (1995). 
[18] E. A. Bartnik, J. A. Tuszynski and D. Sept, Phys. Lett. A 204, 263 (1995). 
[19] Xiao-feng Pang, Phys. Rev. E 62, 6989 (2000); Eur. Phys. J. B 19, 297 (2001). 
[20] D.W. Brown and Z. Ivic, Phys. Rev. B 40, 9876 (1989); 
[21] X. Wang, D. W. Brown and K. Lindenberg, Phys. Rev. B 39, 5366 (1989). 
[22] E. T. Jaynes and F. W. Cummings, Proc. IEE 51, 89 (1963). 
[23] A. J. Leggett et al., Rev. Mod. Phys. 59, 1 (1987). 

[24] W. G. Unruh, Phys. Rev. A 51, 992 (1995); L. Viola and S. Lloyd, Phys. Rev. A 58, 2733 (1998); quantum-ph/980958; 

W. Y. Hwang, H. Lee, D. D. Anh and S.W. Hwang, Phys. Rev. A 62, 062305 (2000); 
[25] G. M. Palma, K.-A. Suominen and A. K. Ekert, Proc. Roy. Soc. Lond. A 452, 567 (1996) [also as |quant-ph/9702001|. 



22 



